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ABSTRACT 

A recent data analysis of the far-infrared (FIR) map of the Galaxy and the Magellanic Clouds 
has shown that there is a tight correlation between two FIR colours: the 60 — 100 /im 
and 100 [im — 140 colours. This FIR colour relation called "main correlation" can be 
interpreted as indicative of a sequence of various interstellar radiation fields with a common 
FIR optical property of grains. In this paper, we constrain the FIR optical properties of grains 
by comparing the calculated FIR colours with the observational main correlation. We show 
that neither of the "standard" grain species (i.e. astronomical silicate and graphite grains) re- 
produces the main correlation. However, if the emissivity index at 100 /im < A < 200 /m is 
changed to ~ 1—1.5 (not ~ 2 as the above two species), the main correlation can be success- 
fully explained. Thus, we propose that the FIR emissivity index is ~ 1—1.5 for the dust in the 
Galaxy and the Magellanic Clouds at 100 fj,ia < A < 200 /im. We also consider the origin 
of the minor correlation called "sub-correlation", which can be used to estimate the Galactic 
star formation rate. 
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1 INTRODUCTION 

Dust grains absorb stellar ultraviolet (UV)-optical light and re- 
process it into far-infrared (FIR), the reby affecting the energet - 
ics of interstellar medium (ISM) (e.g. iHirashita & Ferraral 12002). 
The FIR luminosity is known to be a good indicator of star for- 
mation rate (SFR) in galaxi es jKennicuttl[l99i : llnoue et alJkood : 
llglesias-Paramo et alj [2004). This can be explained if a large part 
of dust grains are heated by UV light from massive stars. The FIR 
spectral energy distribution (SED) of dust grains reflects various 
information on the grains themselves and on the sour ces of grain 
heating (e.g. iTakeuchi etai]l2005t iDopita et all [2005). Dust tem- 
perature, which can be estimated from wavelength at the peak of 
FIR SED, is determined by the intensity of the interstellar radiation 
field (ISRF) and the optical properties of dust (i.e. how it absorbs 
and emits light). Since dust grains absorb UV light efficiently, the 
FIR luminos ity and the dust temperature mostly reflect the UV ra- 
diation field jBuat&XvJI 19961) . 

The FIR SED of the Galaxy (Milky Way) - the "near- 
est" g alaxy - has been investigated by v ario us authors. For ex- 
ample, iDesert. Boulanger. & Pugetl jl990l) and lDwek etail H997t) 
explain the Galactic F IR SED as well as the extinction curve. 
iLi&GreenbergU 19971) provide a unified explanation for the Galac- 
tic polarization and extinction curve by adopting grains composed 
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of silicate cores and organic refractory mantles, small carbonaceous 
particles, and polycyclic aromatic hydrocarbons (PAHs). Although 
there are some variations for assumed species, three grain compo- 
nents are usually assumed: silicate, graphite and PAHs. Besides the 
grain composition and the ISRF, the size distribution also affects 
the FIR SED. The temperature of silicate and graphite grains stays 
at an equilibrium temperature determined by t he radiative equilib- 
rium if the size is roughly larger than 0.01 /im jDraine & Anderson] 
Il985t) . Such large grains are called large grains (LGs). Since the 
equilibrium temperature is usually ~ 15-30 K in galactic envi- 
ronments, the emission peak lies around wavelengths of ~ 100- 
200 fim. If grains are smaller, they are transiently (stochastically) 
heated to a large temper ature and emit photons wi th wavelengths 
shorter than ~ 60 ^tm toraine & Andersoi]|l985l) . Those grains 
are called very small grains (VSGs). The other component, PAH, 
only contributes to wavelengths shorter than ~ 20 /im. Since we 
are interested in wavelengths longer than ~ 60 fim, we do not treat 
PAHs in this paper. 

So far, a major part of works have focused on FIR emission 
from high Galactic-lat itude cirrus or from some limited regi ons in 
the Galactic plane (e.g. lDwek et al.ll997l ; lLagache et aljl998l) . Be- 
cause of the complication in the Galactic plane, which has large in- 
homogeneity with a variety of ISRF intensity and with a large range 
of gas (dust) density, etc., one may expect that a unified under- 
standing of FIR S ED in the Galactic plane is impossible. However, 
iHibi et all (|2006, hereafter H06) show that the FIR colour-colour 
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relation (60-100 /im vs. 140-100 /im) defined by data points with 
Galactic latitudes of \l\ ^ 5° has a clear sequence. A large frac- 
tion (~ 90-95%) of the data points are on a colour-colour rela- 
tion which they call "main correlation". A minor part (~ 5-10%) 
of p oints belong to another sequence called "sub-correlation" (see 
also lSakon et al.l2004l2007h . The sub-correlation tends to be found 
in regions with large radio intensity, which indicates a strong ISRF. 

The "main correlation" of the Galactic plane is also consis- 
te nt with the FIR colou rs of a sample of nearby galaxies compiled 
by iNagata etal] d2002l) (see Fig. 8 of H06). In particular, the FIR 
colour-colour relation (60-100 /im vs. 140-100 /im) of the Magel- 
lanic Clouds is located at the extension of the main correlation of 
the Galactic plane. As stated by H06, this implies that the optical 
properties of grains in FIR are common among the Galaxy and the 
M agellanic Clouds. H 06 also show that the FIR colours predicted 
by Li & Draine d200lh based on "standard" grain optical properties 
( IDraine & Led 19841) cannot reproduce the main correlation. 

Indeed, some observational and experimental results show that 
the FIR grain properties may be different from the "standard" ones. 
Based on the data of the Far-Infrared Absolute Spectrophotome- 
ter (FIRAS ) aboar d the Cosmic Background Explorer (COBE), 
iReach et"al] d 1995b suggest that the Galactic dust has more en- 
hanced emissivity at A > 100 /im than expected before. Some 
ISOP HOT observatio ns have also reached a similar conclusion 
ddeT Burg o et al .120031) . although those observation target dense re- 
gions, where the dust properties may be modified by grain coagu- 
lation to form fluffy aggregate s and/or by coating of ice mantles. 
Recently, Doba shi et al.l J2005I) have shown that the extinction at 
V band derived from the FIR dust emission is higher than that 
derived from the optical Digitized Sky Survey data, suggesting a 
previous underestimation of dust emissivity at FIR wavelengths. 
An e nhanced FIR emiss i vity is also derived f or extragalactic ob- 
jects dDasvra et alj2005T ). lAgladze et al.Ul996i) find experimentally 
that some kin ds of amorphous silica te show emissivities larger than 
those found in IDraine & Leeld 19841) . Moreover, the wavelength de- 
pendence sugge sted by those authors is different from that pre- 
sented bv lDraine &LeeNl984l) . Therefore, it is worth reexamining 
the FIR emissivity of dust by using the new observational results 
provided by H06. 

The aim of this paper is to further investigate the origin of 
the main correlation both observationally and theoretically. Since 
the same main correlation is defined for the Magellanic Clouds, we 
also treat the FIR emission in these galaxies. Furthermore, we ex- 
tend the analysis of H06 to high-Galactic-latitude regions to exam- 
ine if the universality of the main correlation found in the Galactic 
plane holds over the entire Galaxy. Since the ISRF in high Galac- 
tic latitudes is considered to be more uniform than in the Galactic 
plane, the dependence on the Galactic latitude provides us with a 
key to understand the effect of inhomogeneous ISR F. A theoreti- 
cal framework for the FIR SED developed by IDraine & Lil (2001, 
hereafter DL01) is adopted for theoretical analysis of the main cor- 
relation. The sub-correlation is also investigated. 

This paper is organized as follows. First, in section [2] we de- 
scribe the models adopted to calculate the FIR SED. Then, in sec- 
tion^ we briefly review the observational results by H06, present- 
ing also some additional analyses. In section [4] the FIR colour- 
colour relation predicted by our models is compared with the ob- 
servational main correlation. Here, the sub-correlation is also dis- 
cussed. Section|6]is devoted to the conclusion. 



2 FIR SED MODEL 

We adopt the theoretical framework of DL01 to calculate the SED 
of dust emission. Since we are interested in a wavelength range of 
A ~ 60-140 /im, we neglect PAHs, w hich contribute to the emis- 
sion at A < 20 fj,m dDwek et al J 19971 ; DL01). Our assumptions on 
the relevant quantities are described below. 



2.1 Interstellar radiation field 

The ISRF in the solar neighborhood is used as the standard. The 
spectrum of the ISRF in t h e sola r neighborhood is modeled by 
iMathis. Mezger. & Panagial dl983l) . The mean ISRF in the solar 
neighborhood integrat ed for all the solid an gle is denoted as 47rJ®. 
Multiplying it with A, Mat his et al.l dl983l) give the following nu- 
merical fitting in units of erg cm -2 s : 

(A Mm < 0.0912) 
38.57A 4 ,; 4 , 172 



4ttAJ? = 



(0.0912 < A Mm sC 0.110) 
2.045 x 10~ 2 A Mm 

(0.110 < Af.ni «S 0.134) (1) 
7.115 x 10" 4 A^ 6678 

(0.134 < A Mm < 0.246) 
iw\ l , m {W 1 B x (Ti) + W 2 B X (T 2 ) + W 3 Bx(T 3 )] 
(0.246 < A Mm ) , 

where A Mm is the wavelength in units of pm, B\{T) is the 
Planck function, (Ti, T 2 , T 3 ) = (7500 K, 4000 K, 3000 K), and 
(Wi, W 2 , W 3 ) = (1.0 x 10~ 14 , LOxKT 13 , 4.0 x HT 13 ).The 
ISRF intensity J\ is scaled proportionally to J® as follows: 



•A 



xJf, 



(2) 



i.e. the spectral shape of the ISRF is fixed and scaled with \ (x — 1 
corresponds to the ISRF intensity in the solar neighborhood). 

2.2 Optical constants 

We assume a grain to be spherical with a radius of a. The ab- 
sorption cross section of the grain is expressed as 7ra 2 Q a b s (A), 
where Q a bs(A) is called absorption efficiency. Th e optical con- 
stants of silicate and graphite grains are taken from IDraine & Leel 
nm). For the UV regime, we adopt "sm oothed astronomical sil- 
icate" from IWeingartner & Draind d200ll) for silicate grains. For 
graphite grains, two-thirds of the part icles have an dielect ric func- 
tion e — e± and one-third e = en (see lDraine & Leelll984l) . 

2.3 Grain size distribution 

The number density of grains with sizes between a and a + da is 
denoted as rii(a) da, where the subscript i denotes a grain species 
(i — sil or gra stands for silicate and graphite, respectively). We 
assume a power-law form for m(a): 

-K 



rii(a) = da (a min < a < a max ) . 



(3) 



where a m i n and a max are the upper and lower cutoffs of 
grain size, respectively, and d is the normalizing constant. 
IMathis. Rumpl. & Nordsieckl d 19771) find that the extinction curve 
observed in the solar neighborhood is reproduced by a composite 
model of graphite and silicate grains with K = 3.5. We assume 
a = 3.5 A and a m ax = 0-25 fim for b oth graph ites and sili- 
cates dMathis et al.ll977llLi & Drainel200lh . although lLT& Draind 
d200ll) adopt more elaborate functional form. We have confirmed 
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that the difference does not affect the following discussion on the 
FIR colours. 

The mixing ratio between silicate and graphite is uncertain 
and d i fferent even among models for the Gala c tic dust ( Dwek et al. 
1 19971 : iTakagi Vansevicius. & Arimotol 120031 : ILi & Drainell200ll) . 
According to IPeil d 19921) . the ratio between those two species is 
different among the Galaxy, the Large Magellanic Cl oud ( LMC ) 
and the Small Magellanic Cloud (SMC). In the SMC . |Pel ill 9921) 
suggest a silicate-dominated dust composition, while Wel tv et al.l 
d200lh argue against such a composition from a study of the deple- 
tion pattern (but see ILi. Misselt. & Wand [2006). Considering such 
large uncertainty in the grain composition, we present the FIR SED 
of each grain species separately to avoid the uncertainty in the mix- 
ing ratio between the species. The effects of mixing two species are 
briefly mentioned in section |4~3l 

As long as we treat the FIR colours of a single species (section 
12.5b . the constant d cancels out. Thus, it is not necessary to deter- 
mine the absolute abundance of graphite and silicate. The relative 
ratio of d between species is important only when we consider a 
mixture of various species (section [4~3b . In this case, d is deter- 
mined by 



47r 3 r ~K i 

— a Sida da , 



(4) 



where pi is the mass density of species i in the ISM, and S; is the 
material density of sp ecies i. We assume s a g = 3.3 g c m~ 3 and 
s gra = 2.2 g cm -3 jjones. Tielens. & Hollenbachll 19961) . In fact, 
we require only the ratio C s ii/C gra (or p s ii/p gra ) to calculate the 
FIR colors in a mixture of the two species. 

2.4 Heat capacity 

We adopt multi-dimensional Debye models for the heat capacity 
of grai ns. The heat capacity of a graphite grain per unit volume is 
(DL01; iTakeuchi et afll2~ 



2003) 



Cgra — 



T 



2504 K 



y n dy 



-y, f n {x) = —f n (x) 



(5) 



(6) 



Jo exp(y/a;) 

where Nc = 1.14 x 10 cm - 3 dTakeuchi et alj2003h is the num- 
ber of carbon atoms contained in the unit volume of the grain, fca 
is the Boltzmann constant, and T is the grain temperature. The heat 
capacity of a silicate grain per unit volume is (DL01) 

T 



C* si i = (iV a - 2)fc E 



1500 K 



(7) 



where 7V a = 8.5 x 10 cm is the number of at oms contained in 
the unit volume of the grain jTakeuchi et alj2003l) . 



2.5 Calculation of FIR colours 

The total FIR intensity (per frequency per hydrogen nucleus per 
solid angle) of dust emission at a wavelength A can be estimated as 

H(X) = / da— n l (a)7ra 2 Q abs (A) / dT B V (T)—^, 
A min nH Jo dT 

(8) 

where i specifies the grain species and dPi/dT is the temperature 
distribution function of the grains. We assume that the FIR radiation 



is optically thin. Thus, the observed colour directly reflects the ratio 
of the above intensity. Here the colour at wavelengths of Ai and A2 
is defined as I v (\i)/ I v {\2), which is called Ai — A2 colour. 

The distribution function dPi / dT is calculated by the formal- 
ism of DL01 (see their section 4) based on the quantities in sec- 
tions [2j}j2!4] The enthalpy of a grain is divided into discrete bins, 
and the transition probability between the bins is treated by taking 
into account the h eating by interstellar radiation and the cooling by 
emission (see also Guhath akurta & Dr aine 1989). 



3 DATA 

3.1 Outline of the analysis by H06 

The data analyzed by H06 are adopted in this paper. H06 used the 
Zodi-Subtracted Mission Average (ZSMA) taken by the Diffuse In- 
frared Background Experiment (DIRBE) of the COBE. They ana- 
lyzed the data at three FIR wavelengths: A = 60 /im, 100 ^m, and 
140 /im. Then, they examined the relation between two colours, 
J„(60 nm)/Iv(lQ0 Mm) and I„(140 pm)//„(100 fim). Follow- 
ing H06, we adopt the intensity after colour correction for the 100- 
/im and 140-pm bands. For the Galactic plane, they treated regions 
with Galactic latitudes of |6| < 5°. The FIR data of the Magellanic 
Clouds are also analyzed. In order to avoid the uncertainty in the 
subtraction of the zodiacal component, they only select the region 
with 7„(60 /im) > 3 MJy sr _1 . The error of I v is less than 10%. 

In Figure [T^, we show the distribution of the Galactic-plane 
data in the colour-colour diagram. H06 found strong correlations 
between the two colours. More than 90% of the data lie on a strong 
correlation called main correlation, which can be fitted as (H06) 



J* (140 /im) 
J„(100 /im) 



= 0.65 



„(60 /tm) \ 
,(100 /im) J 



(9) 



The main correlation also explains the FIR colours of the LMC 
and the SMC (H06b). There is a nother corr elation sequence, called 
sub-correlation in H06 (see also lHibil2006l) : 



I„(140/im) 
J„(100 /im) 



0.93 



I„(60/xm) V 
/„(100 /tm) J 



(10) 



H06 also examined the dependence on the Galactic longitude 
I, showing that the FIR colours tend to shift downwards along the 
main correlation as / increases from 0° to 180°. Since the ISRF in 
the region toward the Gala ctic centre tends to be higher than that 
in the opposite direction dMathis et all 19831) . the colour shift along 
the main correlation can be caused by the difference in ISRF; that 
is, the FIR colours of dust grains irradiated by a high ISRF tends to 
be situated in an upper part of the main correlation on the colour- 
colour diagram. 

We extend our analysis to high Galactic latitudes (|6| > 5°). 
The procedure of the data analysis and the selection criterion for 
the intensity (i.e. J„(60 /tm) > 3 MJy sr _1 ) are the same as H06. 
Most of the data are taken from |6| < 40°, since the criterion for 
the intensity is hard to be satisfied in higher Galactic latitudes. In 
Figure QJ1, we show the results. Although the data distribution is 
shifted slightly downward compared with the Galactic plane data, 
we clearly see the main correlation. Moreover, the sub-correlation 
is not clear in the high Galactic latitudes. This supports the idea of 
H06 that the sub-correlation is produced by the contamination of 
high-ISRF regions, which tend to reside in the Galactic plane. 
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Figure 1. Far-infrared colour-colour (60 /an — 100 fim colour vs. 140 /an — 100 (im colour) diagram derived from the COBE/DIRBE ZSMA map data, 
(a) The data of the Galactic plane (|6| < 5°), the LMC, and the SMC are shown. The contours show the distribution of the Galactic plane data. The contours 
show the levels where 50%, 80%, 90%, 95% and 99% of the data are contained. The LMC and SMC data are plotted with the filled circles and with the open 
diamonds, respectively. The dashed line and the dotted lines indicate the fittings for the main correlation and the sub-correlation, respectively, (b) Distribution 
of the data in high Galactic latutudes (\b\ > 5°). The levels of the contours and the dashed and dotted lines are the same as those in Panel (a). 



3.2 The main correlation 



3.3 The sub-correlation 



The data presented in Figure[T]indicate that the main correlation is 
robust against the change of Galactic latitudes (6). H06 also show 
that the FIR colours shifts along the main correlation by the varia- 
tion of Galactic longitudes (/). These facts mean that the main cor- 
relation is independent of the complexity of the Galactic structure. 
Since the inhomogeneity of ISRF should be smaller in high Galac- 
tic latitudes than in the Galactic plane, the robustness of the main 
correlation over a variety of b suggests that a detailed modeling of 
galactic structures and of inhomogeneous ISRF are not necessary 
to understand the main correlation. 

As mentioned in the previous subsection, H06 argued that the 
main correlation is produced by a sequence of varying ISRF: the 
FIR colours tend to be located in an upper part of the main cor- 
relation for a high ISRF. However, this was only suggested from 
the longitude dependence but was not theoretically demonstrated. 
Thus, in section [4] we investigate how the FIR colour-colour re- 
lation changes for various ISRF intensities. By examining if the 
observed colour-colour relation is explained or not, we will finally 
be able to test what kind of grain species is consistent with the ob- 
servational FIR colour-colour relation ( section 15711 ■ 

It is worth noting that the main correlation can also explain the 
FIR colours in the Magellanic Clouds. Since the LMC has a face- 
on geometry, the inhomogeneity in the ISRF in a line of sight may 
be small. The fact that the FIR colour-colour relation of the LMC 
shifts upward relative to that of the Galaxy suggests that the loca- 
tion in the main correlation reflects the strength of ISRF, since the 
LMC is generally believed to have a hi gher ISRF than the Galaxy 
(e.g. iFukuiluOOll ; iTumlinson et alj|2002l) . It is surprising that the 
FIR colours of the LMC and the SMC are located on the main cor- 
relation defined by the Galactic dust grains, and it is interesting 
to give a unified understanding of the FIR colours for those three 
galaxies. We also tackle this problem in sections[4]and[5J 



H06 show that the sub-correlation can be reproduced by summing 
two colours in the main correlation. Thus, they argue that the main 
correlation is fundamental. We also follow their argument and con- 
centrate on explaining the main correlation. The sub-correlation is 
investigated in section [52l 



4 RESULTS 

4.1 Standard grain optical constants 

From various observational evidence such as extinction curves and 
interstellar depletion patterns, dust is believed to be mainly com- 
posed of two species: s ilicate and carbonaceous d ust (Mathi sTl990l : 
iLi & Greenber3ll997l: Ijone s 2000; O kada et alj|2006h . For those 
two species, Draine & Led d 19841) proposed astronomical silicate 
and graphite, which reproduce the observational extin ction curve 
and FIR emission spectrum (see also lLi & DrainelfeoOlh . 

Based on the models described in section[2] we examine how 
the FIR colours change by the variation of ISRF (x). Follow- 
ing H06, we focus on the relation between the 60 /im-lOO 
colour [1^(60 ^im)//„(100 /im)] and the 140 /im-100 /im colour 
[I v (140 /im) //„ ( 100 (jm)] . In Figures|2] we plot the colour-colour 
relations of graphite and astronomical silicate for \ = 0.3, 1, 3, 
10, and 30. The discrepancy between the observational main cor- 
relation and the theoretical prediction is significant. Some mecha- 
nism that shifts the colour-colour relation upward and/or rightward 
should be included to explain the observational colour correlations. 

It should be noticed that both silicate and graphite grains trace 
similar tracks on the FIR colour-colour relation although the UV 
absorption coefficients and heat capacities are different between 
those two species. This strongly suggests that the track on the FIR 
colour-colour diagram is determined mainly by the FIR absorption 
coefficients of grains. Indeed, both have dependence roughly de- 
scribed as Qabs oc A -2 . Thus, we suspect that the FIR colour- 
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Figure 2. FIR colour-colour (60 /an — 100 (im vs. 140 fim — 
100 (im) relation for silicate and graphite grains (large crosses and large 
squares.respectively). The different points correspond to x = 0.3, 1, 3, 10, 
and 30 from low to high 140 /im — 100 colours. The same contour as 
shown in Figure [T^ as well as the main and sub-correlations (dashed and 
dotted lines, respectively) are also presented. 



colour relation is not sensitive to various factors other than the 
wavelength dependence of the FIR absorption coefficient. 

The colour-colour relation shifts upwards if the contribution 
from VSGs relative to that of LGs increases, since the emis- 
sion from stochastically heated VSGs contributes significant to the 
emission at A = 60 ji m. Indeed, grains w ith a ~ 50 A contribute 
to the 60 fj,m intensity jDraine & 0120071) . A change of the ISRF 
spectrum and/or of the grain size distribution can change the rela- 
tive contribution from the stochastic VSG emission. First we vary 
the ISRF spectrum. A hard spectrum could relatively enhance the 
emission from VSGs, since the cross section of a V SG is sensitive 
to the change of wavelength in the UV regime (e.g. iDraine & Led 
1 1984) . We calculate the FIR colours by assuming a "hard" spectrum 
without the optical-NIR bump; i.e. we adopt 7.115 x 10 -4 A"™ 6678 
for A Mm > 0.246 instead of the expression in equation {T}. The 
FIR colour-colour relation calculated with this hard spectrum is 
shown in Figure [3^. The colour-colour relation shifts upward and 
approaches the main correlation compared with the relation in Fig- 
ure[2] This is due to relative enhancement of the stochastic heating 
of VSGs, which contributes to the 60 fim emission. However, the 
change of the FIR colours for various \ d° es n °t follow the main 
correlation. Moreover, such a hard spectrum as assumed here is un- 
realistic for high Galactic latitude, where the main correlation is 
also seen. Thus, we conclude that the change of the UV spectrum 
cannot explain the observational FIR colours. 

The enhancement of the stochastic heating occurs also if the 
relative number of VSGs increases. Thus, we expect that the colour- 
colour relation shifts upward if we assume a grain size distribu- 
tion with an enhanced num ber of VSGs. The number of VSGs may 
be increa sed by shattering dJones et al . 1996). The size distribution 
shown bv lJones et al.l dl996) can be roughly fitted by n(a) oc a -4 
(i.e. K = 4.0 in equationfj). We also examine K — 3.75 as an in- 
termediate case. In both cases, we adopt the same values of a m i n 
and a max as assumed in section 12.31 In Figure fjp, we show the 
colour-colour relation for K — 3.75 and 4.0. We observe that the 



60 /im-100 colour is sensitive to the grain size distribution. 
However, the overall trend of the main correlation is not repro- 
duced neither by K = 3.75 nor by K — 4.0. The 60 /xm-100 
/im colour is clearly out of the observed range for K — 4.0. Al- 
though K = 3.75 explains the most concentrated part of the data 
with x ~ 0.3-1, the trend along the main correlation is not repro- 
duced. 

The upper and lower cutoffs of grain size (a ma x and a m in) 
can also be changed. However, the 140 /im-100 /im colour is not 
sensitive to the change of o m i n because it is determined by the 
equilibrium temperature of large grains. It is not largely affected 
by the change of a max as long as a max > 0.1 /im. The 60 /im- 
100 colour monot onically increases if the fraction of VSGs in- 
creases. According to IDraine & Lil d2007t) . grains with a ~ 50 A 
contribute to the 60 /im flux, while a large fraction of the 100 /im 
flux comes from grains with a > 100 A. Thus, as a m i n decreases 
and/or a max decreases, the 60 /im-100 fim colour monotonically 
increases. However, this colour is not sensitive to the grain size dis- 
tribution at a <C 50 A, and the abundance of s uch very small par - 
ticles is constrained better in mid-infrared (e.g. lSakon et ai]|2007h . 
We have confirmed after calculations with our models that resulting 
colour-colour relations predicted by decreased a m i n and/or a max 
are similar to those in Figure[3]and that the observed main correla- 
tion can never be explained by changes of aw n and o max . 

In summary, it is difficult to explain the overall main correla- 
tion by simply enhancing the stochastic heating component, which 
contributes to the 60 fim intensity. Heat capacities and grain optical 
constants are also concerned with the FIR SED. However, as shown 
in AppendixfA] difference in the heat capacity does not change the 
FIR colour-colour relation. There, we also show that different ma- 
terials with different UV-near-infrared absorption coefficients does 
not explain the observational FIR colour-colour relation. We note 
that all the materials examined in AppendixlAlhave the FIR spectral 
index (which is defined as f3 = — d In Q a b s (A)/d In A, i.e. Q a b s is 
approximated to be proportional to A - ' 3 ) nearly equal to 2. We are 
left with changing the FIR optical properties of grains. The follow- 
ing subsection is devoted to examining the possibility that the FIR 
spectral index is different from 2. 



4.2 Dependence on the FIR emissivity index 

A spectral index (5 = 2 is the most frequently used val ue, which is 
supported by some experim ents and observations (e.g. lHildebrandl 
1 19831 : IDraine & Leel Il984l). On the ot her hand, an analysis of 
the FIRAS data by Wright et a I] Jl99ll) indicates that the Galac- 
tic dust has j3 ~ 1.65. Some other observational and experi- 



mental results also indicate that the FIR emissivity inde x may 
be smaller than 2 dReach etalj|l995t lAguirre et alTl2003b P1 The 
range of j3 is roughly from ~ 0.9 to 1.6 as sugge sted for the 
FIR colours of the Galaxy and the Magellanic Clouds ( Rea ch et al.l 
1 19951 lAguirre et alj|2003[). Some spe cies indeed show such a low 
emissivity index dAgladze et al.| [l9961: /3 ~ 1.2 for MgO • 2Si0 2 
(amorphous). Since there is no available information on the opti- 
cal constants of MgO • 2SiC>2 in the wavelengths shorter than FIR, 
we cannot include this species consistently into our models. For- 
tunately, the resulting FIR colours are not sensitive to the detailed 



1 H06 adopt a colour correction with a spectral index of 2 (/3 = 2), but 
even if the colour is corrected by assuming /3 = 0, the difference in the 
colours is less than 10%. Thus, the prediction with j3 = 1 can be compared 
with the results of H06 within an uncertainty of < 10%. 
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Figure 3. Same as Figure[2]but for (a) the hard ISRF and (b) the size distribution of grains with K = 3.75 and 4.0 (lower and upper sequences, respectively). 



wavelength dependence of the UV absorption coefficient as shown 
in Appendix [A] Thus, we adopt the emissivity of astronomical sil- 
icate or graphite at A ^ 100 fim and change the spectral index 
at A > 100 /jm. That is, we choose the following absorption effi- 
ciency for A > 100 fim: 



Q abs (A, 13) = Q abs (A = 100 /im) 



A 



100 jjm 



(11) 



where we adopt the value Q a b s (A = 100 fim) from lDraine & Leel 
dl984h for astronomical silicate and graphite. 

In Figure |4] we show the colour-colour relation for various 
FIR spectral indexes. Here we call the species according to the 
adopted optical properties at A < 100 fim although we modify 
artificially Q a bs at A > 100 fim. We observe from Figure [4] that 
= 1 is roughly consistent with the main correlation if the ISRF 
intensity is between ~ 1 and ~ 10 (/3 = 1.5 may also be permitted 
considering the uncertainty in the data). It is worth noting that the 
colour-colour relation is just along the main correlation. Moreover, 
the region of the largest concentration of the data (see the contours 
in Fig. [4} is reproduced with a standard ISRF % ~ !■ Thus, we 
conclude that a FIR emissivity index of ~ 1—1.5 is successful in 
reproducing the main correlation. 

For the absorption efficiency at A > 100 fim, we als o examine 
a hybrid functional form proposed by Reac h et al] d 1995h : 



Qabs(A) 



(A/Ao) 



[l + (Al/A)6]l/6 



(12) 



which behaves like equation i ll It with j3 = 1 for A <g Ai and j3 = 
2 at A 3> Ai (see also Bian chi. Davies. & Altonl l 1999k We assume 
that Ai = 200 fim, following lReach et alj<199a) . For A < 100 fim, 
we adopt the same absorption efficiency of astronomical silicate 
and graphite as used in section |4~T1 Thus, we set the value of Ao so 
that the continuity of Q a bs at 100 fim is satisfied. In Figure [5] we 
show the colour-colour relation for this absorption efficiency. The 
resulting relation is similar to that in Figure|4]with f3 = 1. Thus, 
the 60 /im-100 fim vs. 140 /im-100 fim colour-colour relation is 
determined by the spectral index around A ~ 100 /im-140 fim 
and is not sensitive to the details of the absorption coefficient at the 
other wavelengths. 



4.3 Silicate-graphite mixture 

In reality, both silicate and graphite contribute to the FIR SED at the 
same time. Thus, it is worth examining the FIR colour-colour rela- 
tion predicted by mixtures of those two species. As stated in section 
12.31 there is a large uncertainty in the mixing ratio among models. 
Thus, we examine a variety of mixing ratios. In Figure|4j;, we show 
the FIR colour-colour relation calculated by a graphite-silicate mix- 
ture. First, the mass ratio of th ose two species i s assumed to be 
Psii/pgra = 1.11 according to Takagi et al. (2003). With this ratio, 
graphite dominates the FIR SED and the FIR colour-colour rela- 
tion becomes almost the same as that of graphite (see Figu re |4p). 
A sim ilar graphite-dominated FIR SED can be seen in lDwek et all 
( 1997). As the fraction of silicate becomes larger, the FIR colour- 
colour relation approaches that shown in Figure |4jt. In order to 
show this, we also examine a silicate-dominated composition as 
proposed for the LMC. We adopt a mass ratio of p s ii/p gra = 7.25. 
The result is shown in Figure |4ji, which becomes relatively simi- 
lar to the silicate colours shown in Figure|4ji. However, it is evident 
that the colours of the mixed population are between those of single 
species. 



5 DISCUSSION 

5.1 Grain optical properties in FIR 

We have found in the previous section that the FIR colours are sen- 
sitive to the FIR absorption coefficient of dust grains and are ro- 
bust to the absorption coefficient in wavelengths other than the FIR. 
Then we have suggested that the FIR emissivity index of ~ 1 fits 
the main correlation. An emissivity index of /3 ~ 1 can also explain 
the submillimetre part of th e Galactic FIR SED without introduc- 
ing a very cold component jReach et al.| [l995)Fl Indeed it would 
be possible to explain the FIR colour by mixing two components 
with different dust temperatures, but the robustness of the main 



2 lLagache et all (1998) argue that the cold component may be due to the 
presence of an isotropic (probably cosmic) background. However, since we 
only select pixels whose 60 fim intensity is larger than 3 MJy sr -1 , the 
contamination of such a background is negligible. 
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Figure 4. Same as Figure[2]but for various FIR emissivity indexes /3 denned at A 100 a™. The filled squares, crosses, and open squares show the FIR 
colour-colour relation for /3 = 1, 1.5, and 2, respectively. In Panels (a) and (b), we assume the silicate and graphite absorption efficiencies at A < 100 fim, 
respectively. In Panel (c) and (d), the colours-colour relation of the graphite-silicate mixture is shown, where the mass ratios between silicate and graphite are 
1.11 : 1 and 7.25 : 1, respectively. 



correlation against the Galactic latitudes would strongly require 
that the mixing ratio between those two components be uniform 
over the entire Galaxy. More natural interpretation is obviously that 
the main correlation reflects the FIR absorption coefficient of dust 
grains itself. Thus, it is important to investigate what kind of the 
FIR absorption coefficients reproduces the observational main cor- 
relation. 

We note that we cannot strongly constrain the spectral index 
at wavelength out of the range treated in this paper and that we do 
not reject the possibility that the spectral index varies as a function 
of wavelength as we examined in equation dl2t . In most of the pre- 
vious works on the Galactic FIR spectra, within the uncertainty, the 
FIR SED at 60 Aim < A < 200 Aim was fitted quite well both by 
/3 = 2 and by f3 — 1 with different dust temperatures. However, 
since the analysis adopted in this paper (H06) only uses the pix- 
els with negligible uncertainty (< 10%), we are able to constrain 
the FIR dust emissivity index more strongly than previous works, 
based on the strong correlation called main correlation. 

Also for extragalactic objects, the FIR spectral index is a mat- 
ter of debate. The FIR-submillimetre (sub-mm) SEDs of dust-rich 



galaxies such as spiral galaxies usually show an excess in the sub- 
mm regime relative to a single component fit to the FIR regime. 
This is usually interpreted as exist ence of very cold dust contribut - 
ing to the sub-mm emission (e.g. IStevens, Amure. & Geadl2~005h . 
However, it is also important to stress that the discussion on the 
amount of very cold dust strongly depends on the assumed FIR- 
sub-mm spectral index. The excess of the sub -mm emissio n could 
be explained by a small spectral index (/3). iKiuchi et al.l ( 2004) 
show that the FIR-sub-mm SEDs of some galaxies can be fitted 
with f3 ~ 1-1.5, similar to our derived values. 

As shown by H06, the main correlation is also consistent with 
the FIR colours of nearby galaxies. This implies that investigating 
the origin of the main correlation is of a general significance to the 
interpretation of the FIR colours in extragalactic objects. 



5.2 Origin of the sub-correlation 

H06 propose that the sub-correlation shown by the dotted line in 
Figure Q] can be explained by mixing two colours in the main cor- 
relation. Such a mixture of FIR colours in a line of sight is called 
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Figure 5. Same as Figure[2]but for the FIR emissivity shown by equation 

M 

overlap effect in H06. H06 also show observationally that the sub- 
correlation tends to be associated with lines of sight with a high 
ISRF intensity indicated by a strong radio continuum intensity. 
Thus, the sub-correlation may be reproduced by a contamination 
of a FIR colour with a high radiation field intensity. 

In order to investigate the overlap effect in our framework, we 
assume a two-component model, in which two FIR colours calcu- 
lated with two different radiation field intensities (y ) are mixed. 
Such a superposition of FIR colours is treated by Ona ka etHEl 
in a more sophisticated manner, but our simple approach 
has an advantage that the physical quantities are easy to interpret. 
We denote the FIR intensity predicted under an ISRF intensity of \ 
as /[/(A; \). Then, we examine a two-component model in which 
the FIR intensity is composed of a mixture of two ISRF \i and X2 
(Xi < X2),I V {\\ {xi, Xa}>/a): 

/,(A; {xi, X2>, fa) = (1 - / a )J„(A; xi) + hh{\ X2) , (13) 

where fa is the contributing fraction of the component with \ — 
X2 in the line of sight. Since the FIR emission is optically thin, 
the contribution from each component is proportional to the dust 
optical depth. Thus, if the total dust optical depth at wavelength A 
is r, far is occupied with a component with x — X2- 

We adopt the FIR absorption efficiency expressed by equation 
i ll It with (3 — 1, since it explains the main correlation (section 
14.21 Figure |4p). We adopt the graphite species for the absorption 
efficiency at A < 100 fim. If we adopt the silicate absorption effi- 
ciency at A < 100 fim, we obtain similar results, and the following 
discussions are valid. We fix xi — 1 as a general ISRF and exam- 
ine various Xz(> 1) an d fa as a contaminating high radiation field. 
In Figure [6] we show the results. We find that the main correla- 
tion shifts rightward (i.e. toward the sub-correlation) on the colour- 
colour diagram if we mix two colours with different x- From Figure 
[6] if about 10% of the total column of a line of sight is filled with a 
high-x region, the sub-correlation is reproduced. In the following, 
we discuss the filling factor assuming that the high radiation field 
originates from young massive stars. 

We consider typical star-forming regions containing OB stars 
as regions with a high radiation field. iHirashita et aU (2001) show 
that a typical Galactic H II region emit ionizing photons at a rate of 



Figure 6. Same as Figure[2]but for the two-component models with %\ = 1 
and X2 = 3, 10, and 30 (lower to upper points). The filling factors of the 
second component (fa) are assumed to be 0.3, 0.1, and 0.03 for the crosses, 
filled squares, and open squares, respectively. 



3.0 x 10 49 s 1 . If we assume that th e mean energy of the ionizing 
photons is 15 eV dlnoue et alJbood) . the luminosity of the ioniz- 
ing photons from the star-forming region becomes 1.9 x 10 Lq. 
Assuming that the luminosity of nonionizin g photons is 1.5 t imes 
as large as that of the ionizing UV photons dlnoue et al . 2000), 
obtain the nonionizing UV luminosity of the star-forming region, 
£\jv, as 2.8 x 10 5 L@, We a ssume that the nonionizing photons 
contribute to the dust heating l lBuat&XuF l996; Inou e et alj200Ch 
and treat the whole stellar cluster as a point source. Then, the UV 
radiation flux (/uv) a t a distance of r from the centre of the star- 
forming region can be estimated as 



/uv 



£vv e" Tuv 
47rr 2 



9.0 x 10" 



10 pc 



— TTJV —1 —2 

e u erg s cm 



(14) 



where ruv is the optical depth of dust grains in the UV regime. 
On the other hand, integrating equation QJ over the UV regime 
(A = 912^1000 A), we obtain the UV flux as / uv = 2.89 x 10~ 3 x 
erg s" 1 cm" 2 (we used 7.115 x 10" 4 A-°- 6678 for A Mm > 0.246). 
Equating this with equation d 14b . we obtain 



X = 31 



10 pc 



(15) 



The UV optical depth (ruv) is related to the UV extinction 
in units of magnitude (Aw) as ruv = -Atjv/1-086. If we take 
A = 2000 A as a representa tive value of the UV wavelengths 
drlirashita. Buat. & Inouel2003h . we obtain Avv = ffa-y (5.52 + 
Rv) = 8.60E B -v for the Galactic extinction curve JPelll 19921) , 
where Eb-v is the colour excess defined between the B and V 
bands (Eb-v = Ab — Ay), and Rv = Av/E B -v (A\ is the 
extinction at a wavelength of A). The ratio between the hydrogen 
column density (both in atomic and molecular form) and E B -v is 
estimated as Nh_ /Eb-v = 5.8 x 10 21 cm" 2 ma g -1 in the Galac- 
tic environment teohlin. Savage. & Drakej fl978). Combining the 
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Figure 7. ISRF normalized to the solar neighborhood value, x, as a function 
of the distance from the center of the star-forming region, r. The solid and 
dotted lines correspond to the results with njj = 10 2 cm -3 and 1 cm -3 , 
respectively. 



above expressions, we finally obtain the relation between ttjv and 
N K asTijy = AW ( 7-32 x 10 20 cm- 2 ). 

iHirashita et alj f200lh assume that the number density of hy- 
drogen nuclei in a star-forming region is typically «h = 10 2 cm -3 . 
Adopting the same number density, we estimate that Nu = nnr 
(with riH = 10 2 cm -3 ) around a star-forming region. In Figure[7] 
we show the ISRF as a function of r. The sharp decline at r > 5 
pc is due to the dust extinction. Although Hirashita et al. (2001) 
focused on H II regions with typical stellar ages less than ~ 10 7 



stars with lifetimes of - 



10° 



r (B-type stars) also emit UV 
radiation (e.g. IHirashita et al]|200lh . Such UV emitting stars with 
long lifetimes, after dispersal of dense regions after ~ 10 7 yr, may 
be found in the field where the ISM density is roughly ~ 1 cm - . 
Thus, we also investigate the radiation field with nn = 1 cm -3 . In 
this case, the radiation field declines as r~ 2 and the dust extinction 
is negligible up to r ~ 100 pc as shown in Figure|7] 

H06 have shown that the FIR colour-colour relation follows 
the sub-correlation in ~ 5-10% of the lines of sight in the Galactic 
plane. The probability (denoted as P(r)) that an arbitrary line of 
sight along the Galactic plane passes through a region with a radius 
of r can be estimated by P(r) — irr 2 sNsF / (~kR\2H c \), where s is 
the path length, is the scale height of the Galactic disk, and Nsf 
is the number of the regions in the entire Galaxy. Taking s ~ Rd, 
we obtain 



P(r) 



0.13 



5pc 



V 10 4 J 



Rd 



10 kpc 



100 pc 



(16) 



This probability is possibly dependent on the line of sight through 
the density structure in the Galactic plane such as spiral arms, ver- 
tical stratification, etc. However, the above probability is valid if 
we take a large number of lines of sight in the Galactic plane as 
H06 have done. The estimate of Nsf in terms of the Galactic SFR 
(or OB star luminosity) is presented in Appendix|B] and the typical 
SFR (or OB star luminosity) of the Galaxy gives Vsf ~8x 10 3 . 



The above probability estimated in equation l !16t is compara- 
ble to the fraction of sightlines with the sub-correlation if r ~ 5 
pc. If the sub-correlation is associated with regions near massive 
stars, the sub-correlation passes through a region whose typical dis- 
tance from young stars is 5 pc. From Figure [7] the radiation field 
at r ~ 5 pc from massive stars is typically \ ~ 10-100. Thus, in 
order to explain the sub-correlation, it would be necessary to con- 
sider a mixture of two FIR colours: one with \ ~ 1 as a normal 
diffuse ISRF and the other with x ~ 10-100 as a radiation field 
near massive stars. This is roughly consistent with the mixture of 
the two components examined above and in Figure|6] 



6 CONCLUSION 

The observational main correlation between the 60 /im-100 /im 
colour and the 140 ^im-100^im colour found by H06 is not ex- 
plained by the previous models. We have examined changes of 
quantities concerning the FIR SED. The most concentrated region 
in the colour-colour diagram and the trend along the main correla- 
tion cannot be explained neither by the change of the ISRF spec- 
trum nor by the the variation of the grain size distribution. The FIR 
colour-colour is not very sensitive to the heat capacity and the de- 
tailed wavelength dependence of the UV absorption coefficient. 

On the other hand, we find that the FIR colour-colour rela- 
tion is sensitive to the FIR absorption coefficient, in particular to 
the emissivity index (5 defined in equation i ll It . The observational 
main correlation is consistent with the theoretical FIR colours if 
we assume that f3 ~ 1—1.5 around 100-200 /im. This is differ- 
ent from the spectral index often assumed for astronomical sili- 
cate and graphite species (J3 ~ 2). Indeed some observational and 
experimental works suggest that the FIR spectr al index is s ignifi- 
cantly smaller than 2 and may be nearly 1 (e.g. Reac h et alii 19951 : 
lAgladzeetaTlll99d) . The difference in the FIR emissivity affects 
the estimate of dust mass and of dust temperature. Thus, we should 
further examine the FIR dust properties not only of the Galaxy but 
also of extragalactic objects by reexa mining FIR data. N ew data 
taken b y the Spitzer Space Telescope I Werner et a l. 2004) and the 
AKARI ( lMurakamj|200l : IShibail2'004l : lMatsuhara et alj2006h . both 
of which have a FIR band with A > 100 /im, are suitable for such 
a purpose. 

We have also considered the origin of the sub-correlation by 
taking into account the overlap effect; that is, the FIR radiation 
comes from regions with a variety of UV radiation field intensity. 
For simplicity, we assume a mixture of two components illuminated 
by different UV radiation fields: one is a normal diffuse component 
with x ~ 1 and the other is an intense field near to young massive 
stars (x ~ 10). As a result, we have found that the sub-correlation 
can be understood by the FIR colour with high x contaminated by 
that with normal x(~ 1)- This explains the observational fact that 
the sub-correlation is associated with the regions with intense ra- 
diation field (H06). Indeed, we have also shown that the observed 
fraction of the data points on the sub-correlation is consistent with 
the total luminosity (or number) of OB stars in the Galaxy. 
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APPENDIX A: MATERIAL DEPENDENCE 
Al Optical constants 

Among silicate species, olivine and pyroxene are often con- 
sidered as promising candidates of astronomical grains (e.g. 
iDorschner et al.l [l995l) . Since the optical constants are different 
among species, we examine the robustness of the predicted FIR 
colour-colour relation against the change of materials. As for 
olivine, the optical constants of two compositions are available: 
Mg 2l ,Fe2-2„Si04 with y = 0.4 and 0.5 for A = 0.20-500 ^m 
jPorschner et alj[l995h . The optical constants at A > 500 /xm 
do not affect our results because of small absorption efficiency, 
but those at A < 0.20 fim may be important, since the absorp- 
tion of UV radiation determines the dust temperature. For both 
y = 0.4 and 0.5, we adopt the optical constants of olivine derived 
bv lHuffman & Stapd ( fl973l) at A < 0.20 ^. 

We calculate the FIR colour-colour (1(60 ^m)//(100 pim) 
vs. J (140 /xm)//(100 /jm)) relation by adopting the above opti- 
cal constants. The other quantities adopted are the same as those 
in the text. As a result, we find that the FIR colour-colour relation 
of the above materials is the same as that of astronomical silicate 
within the error of the colour-colour relation. Thus, olivine is not 
successful in reproducing the main co rrelation. If we adopt th e op- 
tical constants of pyroxene found in Dorschner et al. ( 1995J) (but 
we use the olivine optical constants in bluffman & Stapddl973l) for 
A < 0.2 fim), the main correlation is not reproduced. Thus, as far 
as we examined here, we conclude that silicate does not explain the 
observational FIR colour-colour relation. 

We have examined other materials (FeO, Fe304, Si02, and 
MgO) that can be included in silicate s by using the optical con- 
stants adopted in Hirashita et alj d2005t) . We also investigate amor- 
phous carbon instead of graphite. The optical constants measured 
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Figure Al. FIR colour-colour (60 fim — 100 fim vs. 140 fim — 100 /im) 
relation for various heat capacities. The other quantities adopted for calcu- 
lations are the same as those in Figure [2] The filled squares represent the 
colour-colour relation with the heat capacity given in the text for graphite, 
and the crosses and open squares show the results with heat capacities mul- 
tiplied by 1/2 and 2, respectively. The contours and the dashed and dotted 
lines are the same as those in Figure[2] 



by lEdd dl983h are used (also adopted by Iffirashita et al J l2005h . 
However, all of those optical constants do not change the result- 
ing FIR colour-colour significantly and do not explain the ob- 
servational main correlation. The wav elength depend e nce of the 
absorption efficiency Q a b s is given in Irlirashita et alj J2005h and 
iTakeucMet all (2005) for the above species. 

In our calculation, the optical constants from UV to FIR are 
required. The above materials have such data and are suitable for 
our work. For completeness, we should examine all kinds of pos- 
sible species, but this is obviously impossible. However, the above 
robustness of the predicted FIR colour-colour relation implies that 
the detailed wavelength dependence of the optical constants is not 
reflected in the FIR colour-colour relation. We also note that the 
wavelength dependence in the FIR absorption efficiency (emissiv- 
ity) is roughly proportional to A~ 2 (i.e. j3 ~ 2) for all the above 
species. As shown in text, the FIR colour-colour relation is sensi- 
tive to the FIR emissivity index. A similar wavelength dependence 
of FIR emissivity is the reason why the above species produce sim- 
ilar FIR colour-colour relations. 



smaller heat capacity, since the temperature of a VSG rises more 
when a photon is absorbed. At the same time, however, the grains 
rapidly cool if the heat capacity is small. Those two effects com- 
pensate, and the 60 fim-\QQ fim colour does not drastically change. 
Moreover, since the contribution from LGs in radiative equilibrium 
at the 60 fim band becomes larger for higher ISRF, the difference 
in the 60 fim-100 fim for \ ^ 3 is negligible. Thus, we conclude 
that the change of grain heat capacity does not affect our results in 
the text. 



APPENDIX B: THE GALACTIC STAR FORMATION 
RATE 



In section [5T2l we introduce the number of regions with young mas- 
sive stars, Nsf, to explain the fraction of sightlines associated with 
the sub-correlation. The number of such regions reflects the recent 
star formation activity of the Galaxy. Thus, we relate A^sf with the 
Galactic star formation rate (SFR). 

We assume that the Galactic disk has a radius of Rd with 
thickness of 2Hd (we assume Rd = 10 kpc and Hd = 100 pc). 
For simplicity, we assume that the regions with massive stars dis- 
tribute uniformly in the Galactic plane. Then, Nsf can be estimated 
as N SF /(nR 2 d H A ). 

The total UV luminosity in the entire galaxy can be estimated 
as A^sf^uv, where £\jy is the UV luminosity of a region with ma s- 
sive stars. This can be converted to the SFR as dlnoue et al.ll2 00oFI 

A^sf^uv 



SFR = 



1.8 x 10 9 L @ 



■ M, 



yr 



(Bl) 



Since Itjv ~ 2.8 x 10 5 L , we obtain SFR = 1.56 x 
10' 4 Ar SF M Q yr" 1 . 

The total lum inosity of the Galac tic OB stars is estimated to be 
~ 2.3 x 10 9 L© dMathis et alJll983h . Equating this with ATsf-^sf, 
we obtain JV SK ~ 8.2 x 10 3 (or SFR = 1.3 M yr -1 ). If this 
number is used in equation l l 1 6t . we obtain the probability that a 
line of sight is associated with a high radiation field of \ ~ 10. 



A2 Heat capacity 

The difference in the heat capacity may affect the temperature of 
the VSGs. The LG temperature is not affected by the heat capacity 
since it is determined b y the equilibrium tem perature independent 
of the heat capacity (e.g. lDraine & Lejl984h . Thus, we expect that 
the difference in heat capacity may cause a variation of the 60 
intensity. Here we test heat capacities whose values are 1/2 and 
2 times of the value adopted in the text for graphite. The resulting 
colour-colour relations are shown in Figure lATl As expected above, 
the 140 ^im-100 fim colour, which reflects the equilibrium temper- 
ature, is not affected by the change of heat capacity. This also sup- 
ports our assumption that the contribution from VSGs is negligible 
at A > 100 fim. The 60 /im-100 flux ratio becomes higher for 



3 In deriving equation IB 1 h we identified L non i on = 0.6lJ^ in Inoue et 
al. (2000) with AT S f4jv- 



